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ABSTRACT 


A ooi^t«r siaulatioa has bssn dsvslopsd to asssss passangar survival 
during tha post<-crash avacuation of a transport catagory aircraft vhan fira is ' 
a aajor thraat. Tha coaputar coda, FTREVAC, co^putas individual passangar 
axit paths and tiaas to axit, taking into account dalays and congastion causad 
by tha intaraction aaong tha passangars and changing cabin conditicxis. Si^la 
Bodals for tha physiological effacts of tiia toxic cabin ataosphara are 
included with provisitm for including aora sophisticated aodals as they ba c oa s 
available. Both wida-body and standard-body aircraft nay be siaulatad. 
passenger characteristics are assigned stochastically froa eiqieriaentally 
derived distributions. Results of siaulations of evacuation trials and 
hypothetical evacuations under fire conditions are presented. 
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SBCnOM t 
nniKXxicnoii 

Oni^rsity of Oayton R*s«ardt Institat* (ODRI) sdutiats hava daaalopad 
a coaputer aodol, PIREVAC, to sianlate passangar awcoation froa a ganaric 
aircraft, %fith provision for a post'-crash scanario ineloding fira* 

1.1 H(X)EL OBJECTIVE 

The aodel objective is support for the study of: 

1 ) The effects of fire-induced toxicants on tiae required for passenger 
evacuation and, hence, probability of passenger survival; 

2) the effects of aircraft design and aaterials on tias required for 
passenger evacuation (both with and without fire); and 

3) evacuation procedures. 

1.2 MODEL CStGAMIZATICXI 

The aodel has three logical aodules. they are: 

1 } C^bin Environs«nt Nodule (CEM) - The CEM describes a two- 
diaeps-onal cabin environaent as a function of tiae. 

Environaental factors include cabin configuration (placeaent of 
seats, aisles, doors) and the effects of fire (teaperature and 
concentrations of toxic gases). Possible future additions to 
this nodule would provide for a three-dinensional environaent 
(the addition of height to the present length and width), as wall 
as the inclusion of sao)ce and crash debris. 

The present coaputer code iaplenents the CEM by reading data 
files whid) provide the tiae-dependent cabin description. These 
data files can be either derived froa test data or the output of 
nathenatical node Is of fire such as MacArthur's (UDRX) DACFIR 
aodel developed for the FAA (see Reference 1). 


1 



ORIQINAl. PAGE IS 
OF POOR 0^>ALrrY 

2) Ttim Huaan Factor Hodula (RFN) - Thm RFM caleulatas tha physiolo- 
gical offsets of firs-rslatod toxicants on huaan as capo bo- 

ha viors* 

tho approach takon is tho oso of a tnuMn rosponoo factor 
which is uaod to oodify passongor oowsosnt paraastors and which 
varios with tho history of cabin conditions. Hiis rosponso fac- 
tor is carrontly calcalatod using tho coneopt of a "Fractional 
Incaoacitation Ooso", following tho Idoas of Sarkos and Crano 
(soo Rsforoncos 2 and 3). tho fractional incapacitation doso 
roorosonts a passongor *s accunulation of toxicants (hoot* gasos) 
as a fraction of tho doso roquirod to incapacitato that 
passongor. Ef facts of tho sopors to toxicants considarod aro 
aapuawd additiws. 

3) Passongor Cgross Nodulo (PEM) - tho PEM sisulatos passongor noso- 
aont. Each Passongor is assi^od an "optinal" (soo soction 4.2) 
oxit routs through tho aircraft configuration and proooods along 
this routs sxibjsct to interaction with other passengers and 
asA>ient conditions. Exit routes can be updated to reflect tho 
chanoes in tho cabin onwironoent, as rolayod fro* tho CEM. 
Passonqor ao^oent behavior will change as dictated by the HFK. 
Passenger position is displayed by "snapshots”, graphical output 
representing the aircraft interior as a function of tine. 

Figure i shows the interaction between the three logical nodules in 
terns of the data flow between then. The nodel is a clock-driven sinulation; 
that is, at given tine increswnts the cabin environsont, passenger physical 
condition, and passenger positions are updated. Mote that tho update incre- 
nents need not be the sane for all phenonena, e.g. , passenger positions can be 
calculated nore often than the cabin atnosphere is updated. (See Appendix A, 
card type A for detailed definition of update tine paraneters.) 

The next three sections of chis report will exanine the three logical 
modules in detail. 
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Figure 1 


Data Flow Between Logical Modules of FIREVAC. 
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CABIN EHVntONMEtlT MODOLE 
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As can ba s««n froa Pigara i, tbm CEM inrocassaa two priaary data atme- 
turas, tha cabin confi^ration and tha cabin ataoaphara profila* nia saction 
will <tescriba tha coapatar iaplaaantation of tha CEH in taras of tha braakdom 
of tha priaary data structuras into thair ooaponant parts and tha procassing 
functions which opera ta on those aoiqwnants* 

2.1 THE CABIN C0NPI60RATI(»1 

Tha cabin oonfigaration aodal inputs describe tha urcraft aa a set of 
nodes or boxes* These nodes can represent saatSr aisle space, exit doors, 
aircraft skin or exit slides. Each node is assuaed large enou^ to hold a 
single passenger. Node location is defined %rith a two diaensional row<-ooluan 
coordinate systea. Typical aircraft configurations are shown in Figures 2 and 
3. The nuabered nodes represent seats with passengers, e.g. , row 2, coluan 3 
is the seating (initial placeaent) assignaent for passenger 1 in both 
Figures. 

Nodes are classified according to the type of space they represent 
(i.e., seat node, aisle node, etc). Passenger speed is then defined in teras 
of the time required to traverse a given node type. 

For each distinct node type the aodel requires the mean and standard 
deviation, the maxiaua and the adniaon tiae of passenger aoveaent through that 
node type. Individual passenger node aovesMnt times are stochastically 
assigned by assuaing norxtal distributions within aaxiaua and miniaua (see 
section 3). The aodel has been exercised using time data supplied by 
J. Gillespie of the FAA as described in Reference 4. The test case was a X-9 
evacuation, the data based on a test deaonstration with 144 subjects conducted 
August 2, 1975. 

The values for interior node types (1-3) were not obtained froa the DC-9 
evacuation. To quote froa the above referenced Gillespie report (p.7) 
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Figure 2. Sample Aircraft Configuration 
with Passenger Placement. 
(B767-200) 
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Figure 3. Sample Aircraft Configuration 
with Passenger Placement, 
(DC-9, Series 180) 
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Data on passonnr aowiMnt within tha aircraft ara 
basad upon tasts dona at CAMI with slxtaan tast 
subjacts ona at a tlna as thay rosa froa thalr saat, 
procaadad to tha alalar and aovad d own tha alsla to 
tha axit* Tharafora, this data doas not includa any 
Intarfaranca af facta batwaan pasaangars as thay sows 
into tha alslas* It has baan notad that bottlanacking 
occurs at tha axlt doors In tast casas run with tha 
data* It was laposslbla to obtain passanqar aowanant data 
Inslda tha aircraft from available evacuation fllns. 

Ralavant data in seconds is as follows: 


NOTE TYPE 

_X 

S 

MAX. 

MIM. 

NODE 

DESCRIPTION 

1 

.253 

.03 

.3 

.2 

row aisles 

2 

.253 

.03 

.3 

.2 

column aisles 

3 

.933 

.106 

1.3 

.75 

seat 


4 

.96 

.33 

2.0 

.5 

Type 

I exit doorway 

5 

2.33 

1 .42 

7 . 

.7 

Type 

III overwing doorway 
(smooth flow) 

6 

1.54 

.0 

2.5 

.6 

Type 

I exit slide 

7 

2.97 

.0 

4.2 

1.1 

Type 

III overwing slide 
(smooth flow) 

8 

1.72 

.87 

5.6 

.7 

Type 

III overwing exit door 
(erratic flow) 

9 

2.09 

.0 

4.0 

1 .4 

Type 

III overwing exit slide 
(erratic flow) 


Where x * Mean tine through node (sec). 

s > Standard (teviation of time-through-node distribution (sec). 

Max. & Min* * Maximum and minimum times through node (sec) (i.e., we 
use a clipped Gaussian distribution. ) 

Gillespie maices a distinction between smooth and erratic flow for the 
ovarwing exits because passengers had to be forced by crew members to maintain 
reasonable traffic flow to one of the overwing exits. It should be noted that 
in our data sets S « 0 for the exit slides (nodes type 6, 1 , 9) whereas the 
Gillespie data includes the experimental data for standard deviations. Our 
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rationale la the assunptlon that passenger behavior on the exit slides Is more 
a function of gravity than of those passenger physical characteristics which 
determine behavior interior to the cabin. 

Also required as input Is the mean time to open a door. The data used 

were: 

Type I exit ■ 10.4 secs. 

Type III ovsrwlng exit - saiooth flow - 12.8 secs. 

Type III ovarwing exit'*erratlc flow » 15.8 secs. 

The model has two variables to describe exit status: IXITPP and IXITRL. 
The first of these defines passenger perception of the exit status, l.e., 
whether the passenger believes the exit will be open when the passenger 
reaches it and, hence, is a good exit target. The second provides the physi- 
cal reality of exit status, l.e., Aether a passenger can actually egress 
throu^ a given exit. IXITPP Is used In the determination of passenger exit 
path and IXITPL is used In the simulation of passenger movement (see section 
4). 


The input variables which ooo^rlse the cabin configuration are listed on 
the B and C card types in Appendix A. 

2.2 THE CABIN ATMOSPHERE PROFILE 

The cabin atmosphere is defined as a set of toxicant values 
(temperature, toxic gas concentrations) which are a function of both time and 
cabin position. Variation in time is achieved by updating the toxicant values 
at set intervals (as defined by an input parameter). Variation as a function 
of cabin position is achieved by assuming each cabin node has its ovm 
atmosphere, that is, a complete set of toxicant values is assumed for each 
cabin node. 

Toxicant values are derived from an input file which represents a 
series of sample readings of those values at various cabin locations. This 
input file can contain experimental data or the output from a mathematical 
model of fire. In either case, the input file may not provide the complete 
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cabin atnoaphara; sampla reading times may not coincide with model atmosphere 
update times, end there is no guarantee that each cabin node will represent a 
sampling location for each (or any) toxicant. 

Ihis means the model must expand on the data provided by the atmosphere 
input file and approximate whatever values may be missing. Each model 
atmosphere update performs a two-step approximation process. The first step 
provides temporal variation, the second spatial variation. 'Riis two-step pro- 
cess mirrors the format of the atmosphere input file. Each input file is 
organized into data sets. Each data set contains one reading for each toxi- 
cant at each of that toxicant's sampling locations. Each reading in a given 
data set is assumed to have been taken at the same time, as measured from 
clock time « 0 for the model scenario. (Note that clock time « 0 is not 
necessarily the start of passenger movement. See Section 4 for details of 
passenger movement) . 

Data sets are aissumed to be in chronological order. At any time during 
the simulation (in particular at atmosphere update times) the model is pre- 
sumed to have read (and stored in memory) two data sets; the first repre- 
senting a time less than or equal to the simulation clock time, and the second 
greater than that clock time. Variation in time for the sampling point values 
is achieved by linear interpolation between appropriate values of the two data 
sets. This is the first step of the cabin atmosphere approximation process. 

It provides a temporally complete profile of toxicant values at toxicant 
sampling locations. 

The next step in the cabin atmosphere approximation process is tc find 
values for each toxicant at each cabin node for the given simulation time. We 
use a weighted average of some (or all) the sampling points. This method was 
chosen for ease of implementation and for the generality it offers. It 
requires no assumptions as to the locations or r-'.-nber of sampling points. 

For a given toxicant: 

Let Px i » 1 , n represent the sampling points 

Let P represent a point of interest 
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Ii«t r«pr«s«nt the toxicant mloa (taaparatura, gas con- 
centration, ate.) at Ttw problaa is to calenlata 

a aalua, v, for the point P. Let raprasant tha 
dlstanca froa P to d., (i.a., if P is cabin noda 
(n ,B^) and P. ■ (n_ , a^ ) than d, • f(n^ - n_ 

(. .2 i'>i . ‘ i 

) •> 

i 

If d^ > 0 for any i than P is a saapla point and no approxiaation is 
raqoired. Otharwisa, lat: 


V 


I 


d.<R 

1 



V. 

X 


where: 


R definas the region of approxiaation (i.e., if d^ >R then 

V. does not contribute to V) 

1 


and a 


^ (d ) 
d.<R i 
1 


Hote that as the distance froa P to P^ (tecreases or increases, tlw value 
nakes a correspondingly lesser or greater oontribution to the weighted 
average v. 


The function, Vd^^, used above is not the only choice of weighting fac- 
tor, and was selected for siaplicity in the absence of any criteria to favor 

another choice. Should future experience reconaend another function (as for 
2 

instance l/d. ) the aodel could easily be changed. 


Note in the above if there is only one saapling point, P^, vrithin the 
region of approxiaation, then V « v^. in this instance we have aerely taken 
the value of the closest point. If P is on a line between two points P^ and 
P^ and those are the only two points within the region of approximation, then 
the abowe is equivalent to linear interpolation. 


Ne currently have two sets of cabin atsosphere data. The first of these 
is derived froa FkA C-133 Fire tests as described in reference (see 
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Raf«r«nc« 2) (sM Pig. 4). w* haw usad this data to provlda an ataoaphara 
for vida-body slanilatlons (soaelficaily In a Boalnq 767). UnfortunataIy» ^Is 
data la glwn for only a singla saapUng location. Mora ooaplata data la 
a wi labia froa a seriaa of NASA flra taata (saa Rafaranca S). Thla data aaa 
uaad to provida aultipla saapllng ooint data for both taa^ratoraa and gaaaa. 
Fig. SA shows tha tast oonfiqnratlon with data oollaetion points and Fig. SB 
showj tha oorrasponding aodal saapling point locations for a DC-9 cabin oon- 
f igaration. 
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FAA C-133 Fire Tests (Sarkos, 1982) 


I 



Figure 4(a). Hazards in Aft Cabin Produced by 
Burning Interior Materials. 
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(b) FIBEVAC Representation of 
Test Cabin. 

Figure 5. Continued. 
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SECTION 3 

THE ROMAN FACTOR MODULE . 

As fiqars i shows, sn HFM input is ths toxicant data of tits CEM, and ths 
HFM output is tiSM rsquirsd for passsn^sr noussMnt, as dsgradsd by ths offsets ^ 
of thoss toxicants. Ths CEM producss an ataosphsrs for sach nods in ths 
aircraft. At giwn tine intsr^ls (as detsminsd fay an input paraiMtsr) ths 
HFM sxaninss each passenger, updating his huaan response factor with respect 
to the toxicants present in the no<te he currently occupies. That hua»n 
response factor is then used to nodify his speed by altering both the tine 
required to trasel through all node types and his reaction tins - the tine 
required to notice a target node is '«ncant. (See section 4 for the details of 
node-to-node noesawnt}. 


3.1 THE FRACTIONAL INCAPACITATKNI IX)SE 

The hunan response factor, R^, is related to F^, the concept of a frac- 
tional incapacitation dose (see References 2 and 3) by: 


As explained by Sartos et. al. (Reference 2), the F^ concept is a hypothetical 
hunan survival eodel whose purpose is to assess the relative isportance of 
each cabin fire hazard: 

The survival nodel described ... is hypothetical. Its sain 
purpose is to provide a mans of predicting the tim-of- 
incaoaci ration within a fire enclosure, based on masurenents 
of elevated tenperature and toxic gases concentrations which 
change, in som cases substantially, with tine. Thus, it is a 
tool for reducing a fairly large nunber of somwhat abstract 
masuremnts into a single, cogent paramter: tim-of- 

incapacitation or the hypothetical tim at which an individual 
can no longer escape fron a fire environnent. How well the 
nodel relates to actual escape potential is unknown ud, 
realistically, cannot be determined. It is known that segmnts 
of the nodel are deficient for lade of available information. 

For example, no data exists on the effect of irritant gases 
(e.g., HCL, HF) on acute human escape potential. (FAA has 
sponsored new research at Southwest Research Institute to 
detemine "the threshold concentration for escape impaiment by 
irritant gases (HCL and acrolein, initially) using a nonhuman 
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priaate aodel amd a relevant behavioral task that can be 
extrapolated to nan.") Thus, the HCL and HF incapacitation 
doses utilized in the aodel are siaply based upon extrapolation 
froa threshold liait values (TLV's) for an 8-hour work 
environaent. Confidence in the aodel is greater for the prediction 
of the relative escape tiae between tests on different aaterial 
systeas tihan on the prediction of absolute escape tiaes. 

(Proa Reference 2, pp. 6-7) 


The (diief virtues of the Fq concept froa the point of view of our aodel 


are: 


1 ) The fact that it does reduce a large nuaber of abstract 
aeasureaents into a single paraaeter, and, hence, one that 
can be easily ^plied to passenger behaviors; and 

2) The fact that it allows for the cuaulative effects of the 
atsmsphere, thus aullowing the passenger's sho?t-tera toxicant 
exposure history to affect his probability of survival. 

The UDRI iaplemntation of the Fq concept makes the assumption that 
yields not only the tine to incapacitation, but also a measure of partial 
impairment, e.g., Fq « 0 no impairment, Fq > 1 »> complete incapacitation, 
Fq 3 .5 >> 50% incapacitation, i.e., passenger speed is decreased by a factor 

of 2. 


The present computer im^'>lementation defines F^ as: 


Fo(t) 


N 

\ 

L 

n»1 


3.61 


M 

^ I 

i 


c. 

i,n 


At 


where: 

Fj^(t) > the Fractional incapacitation Dose accumulated at time t 

At = the time increment {in minutes) for the n^ interval 
(not necessarily the same for all n) 

N = the number of time increments to time, t 

N 

i.e., t » V At 
n-1 " 

T^ = the temperature (*C) at time t 
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a 

* 4.1 X 10 statistically deriv«d proportionality constant (so* 
Rsferance 2) 

0^ ■ th« incapacitation dose of the i^ constituent (ppa^see) 

th 

C. « the concentration at the i constituent (ppn) 
i,n 

The constituents currently under considetction are: 

CO » 1.44 X 10® (ppB«sec) 

HF D. « 6.84 X 10^ 

^ 5 

HCL D. « 1.44 X 10 

^ 4 

HCN D. - 2.88 X 10 

1 

Note that the aboue ev|uation assuaes that all effects are additi^. If an 
individual could simultaneously absorb the incapacitiation dose of tvo dif» 
ferent toxicants the equation would give him an equal to 2; however, the 
computer implementation of the equation imposes an upper limit of i on any 
individual's F^. 

Note, furthermore, that this form of F^ does not taUce into account 
individual passengers' respiration rates or body masses. Also, there is no 
consideration of oxygen deprivation or the physiological and/or psychological 
effects of smoke. Future model enhancement should provide for further 
deterioration of passencer speed due to the blinding effects of dense smoke. 

3.2 MODIFICATION OF PASSENtJER BEHAVIOR 

At present, passenger speed (defined as time to move through a node) is 

initially assigned stochastically, using a set of random Gaussian deviates. 

For each passenger P, his speed of movement is determined by first randomly 

selecting a value Z from a standard normal distribution (mean 0, standard 

P « 

deviation 1 ) . Each node type, n, has its associated and mean and 

standard deviation or the time required to move through that node (as defined 
by the CEM). Then t^,p (the time required for the pth passenger to move 
through the n^jj node) is initially; 

t » X + S , if tmin < X + Z S < tmax 
n,p n Pn n— n pn— n 
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and if: 

^ + ZpSn < t^ttn <^,p * 

X„ ♦ 2pS„ > t»ax„ than t„^p - taax„ 
where: 

tainjj is tiie ainiauB tiae allowed for aoveaent through a type n 
node and taaxn is the aaxiaua such tiae. This restriction on ^e 
range of t„,p is required to avoid the aberrations which could 
arise froa blindly fitting a continuous noraal distribution to 
experiaental data (e.g. , If « 2.33, » 1.42, then “ >2 

would give t > .51 without a ainiaua nmge restriction, 
whereas, tain„ > .7 results in tj^,p « .7 and avoids negative tiae 
of aoveaent) . 

The initial values for tn,p are assigned under the assuaption tiiat Fq « 0 and, 
hence, Rg » i . At each HFM update interval, ^e values for t^,p and each 
passenger's reaction tiae are divided by Rg, Hence, when Pd » 0»tn^p is 
unchanged; when Fo * 1/ t,i,p i-s infinite and the passenger is unable to aove. 
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SECTION 4 

PASSENGER EGRESS NCM^EZ. 


The PQI coopletes the cycle of nodule interaction illustrated by 
Figure l « It accepts data on cabin conditions froa the CEM euid data on 
passenger behavior from the HFN, uses that data to simulate passenger nove- 
ment, and returns passenger position data to both the CEM and the HFM as well 
as producing the ocdel's graphical and summary outputs* 


The simulation of passenger movement assumes that at any given time each 
passenger has a Icnown exit path, a sequence of nodes beginning with the 
passenger's current position and terminating with an exit. Each passenger is 
examined each update of the PEM to determine whether that passenger satisfies 
the criteria for movement into the next node of his exit path. 


4.1 PASSENGER MOVEMENT CRITERIA 


For each PEM update for each passenger, p, the simulation logic requires 
variables: 


T * current time according to simulation running clock 


N = row, column location of node currently occupied by p 
in 

N =• row, column Location of p's target node (next node in p's exit 
to 

path) 


T “ time N was vacated (i.e., simulation clock time at which the 

empty to 

last passenger to occupy N left) 


T - time recruired for p to move though N (sec'c) 
in in 

T - time reguired for p to move through N (sec's) 
to to 

T * time of o's last move (clock time) 

lastmove 

REACT * time reguired bv p to notice that N . , , , 

to is empty ( sec ' s ) 


STATUS ■ the passenger number of the current occupant of node n, 
n 

e.g. , STATUS^ » p. STATUS » 0 if node n is empty, 
in 

i.e., D will be unable to nove unless STATUS = 0. 

N 

to 
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Pass«ng«r p is considered to have had enough time to move when: 

1) T — ‘^lastaove Zm (i/2) and 

2) T - Teapty >, REACT 

This technique assuams movement is from the center of the node currently 
occupied to the center of the target node. Condition (2) provides a simula- 
tion of reaction time delay. 

As a result, the model shows passengers in a tightly packed queue moving 
in a shuffling fashion, where movement is jerky and the movement of a given 
passenger is dependent on that of passengers in front of him. Passengers in 
less crovfded quarters are modeled as accelerating to their speed of movement 
and maintuning that speed. This is because in that case, target trades have 
few passengers in them as blockers and, hence, condition (2) is met virtually 
every time condition (1) is met. 

Both the time reqiiired to move through nodes and REACT are modified in 
the HFM by the human response factor. At present, REACT is an input parameter 
and a single value is assumed initially for all passengers. Given experimen- 
tal data, this single value could be replaced by a mean, standard deviation, 
maximum and minimum as with passenger node times. The REACT parameter also 
requires further study to determine proper values for the simulation of panic 
situations, in which passengers would probably be pushing and shoving, amd 
hence, packed nore densely in their exit queues than >rauld be the case in am 
orderly evacuation. 

When passenger p has .am empty target node * Nto, and is found to have 
enough time to move, the nodes adjacent to are examined for other conten- 
ders, other passengers who also meet the above given criteria for movement 
into at the current simulation clock time. If auny other contenders are 
found, priority is given to the passenger who has been waiting longest for the 
given target node. This procedure could easily be replaced in the computer 
implementation if another priority scheme is found to yield a more realistic 
simulation. Other possible decision procedures which have been considered 
are: 
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1 ) Sp««d *• fastest passenger has ri^t’*of-way; 

2) Size - biggest passenger has ri<^t-of-*way; 

3) "Chivalry" - male passengers allow female passengers right-of-way; 

4} "Parental Agressiveness" - passengers identified as carrying small 

children have right-of-way; and 

5) Random draw. 

In terms of the simulation objections, the priority scheme used is not 
as imoortant as such model parameters as passenger speed or REACT, because 
the priority scheme has more effect on which passengers escape than it does on 
how many escape. 

4.2 PASSENGER EXIT PATHS 

The assumption is made that at any time in the simulation, each 
passenger is following a set path to egress, rather than looking ahead only 
one move at a time. In order to have the passenger's movement respond to 
cabin conditions, these set exit paths are updated periodically. In order for 
such updates to make sense, the choice-of-path is dynamic; it reflects 
changing conditions as reported by the CIM. 

The model's exit oath algorithm allows for determination of an "optimal" 
route from a passenger's present node position to the closest exit perceived 
as open by that passenger. The path is optimal in the sense that the 
algorithm calculates a "distance" from the passenger to all possible exits and 
chooses the closest one (as defined by chat "distance"). The algorithm's 
flexibility lies in the number of wavs in which it is possible to measure the 
desirability of a path; at present, the time of movement from one node to the 
next; and the difficulty (or impossibility) of moving through blocked nodes 
is considered. 

The exit route is selected by viewing the aircraft as a digraph 
(directed graph). The node centers are graph vertices and the paths from one 
node to the next are viewed as edges. Finding the exit path is, thus, the 
problem of finding the shortest oath from a specified vertex (passenger's pre- 
sent position) to a'lother specified vertex (open exit). 
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The algorithm used is due to Oijkstra (see for example Reference 6) and 
makes use of the length, <^i,j distance or weight) of the directed edge 
from vertex i to vertex j. This length or distance will determine the desira- 
bility of moving from one node to another and can be defined in a number of 
ways. The only restrictions, the algorithm places on the definition of 
dj^^j are; 

<^i,j2.0 

di i ■ 0 Vi 

dj^^j » • if there is no edge (or path) from i to j 

The current model implementation defines the metric <^i,j for each 
passenger, p, in terms of that passenger's node movement times. The presence 
of other passengers in nodes along a potential exit path is considered a 
possible impediment and affects <^i,j by adding a term designed to represent 
the delay created by waiting for those passengers to move. 

The form used is; 

di,j « (1/2) Ti,p + (1/2) Tj,p > Tj,b 

where: 

Tj^^p » time required for passenger p to move through node i 

T j ^ p » time required for passenger p to move through node j 

0 if there is no passenger in node j 

- 

time required for the blocking passenger (the passenger in 
node j ) to move through node j . 

A sample calculation is shown in Figure 6. In both cases a and b 
passenger 1 is determining his closest exit. Nodes A1 - A6 represent aisle 
nodes, through which passenger 1 can move in .25 sec.; nodes S') - S6 represent 
seat nodes, through which the passenger can move in .9 sec., and nodes E1 and 

E2 represent exits (of the same type), with 1.0 sec. as required time of move- 

ment. In case a, passenger l is assumed alone in the portion of the aircraft 
represented. The digraph representation shown is labeled with the values for 
d between each at the nodes represented. As shown, the closest exit to 
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Aircraft configuration 
with passanqar 1 in AS 
whosa novanant tiiMS ara: 


Digraph rapraaantation 


Bara i ' s axlt path to 
S1 AS A2 A1 SI haa 
through S. (aaat noda) ■ .9 sac. iistanca 1.375 

A. (aiala noda) • .25 sac 

E. (arit noda) « 1.0 sac Exit path to S2 

AS A2 A3 A4 E2 haa 
diatanca 1.375 



If a paaaangar, 2 , is addad in noda At with novaaant tisMS 
through S. « 1.0 sac 

A. • .3 sac 

E. ■ 1.5 sac. tiian tha apparent distanca (froa I's 

parspactiva) froa A2 to A1 
changaa aa shown balow. 
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Hanca 1 ' s exit path to El 
now haa distanca 1.42S 
whareas tha distanca to S2 
is unchanged. 



(b) Two Passenger Casa 



LEGEND ; • passenger 1 

passenger 2 

Figure 6. Illustration of Exit Path Distance Calculations. 
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p«ss«n 9 «r 1 Is El and his sxlt path is A5A2A1E1. In cass b» a sscond, slowsr 
passsngsr is addsd in nods hi . this changss passsngsr l * s psrcsption of ths 
distancs froa A2 to hi by .3 ssconds ths tias rsquirsd for 2 to travsl 
through nods hi . Now passsngsr 1 bslisvss E2 is ths clossst sxit, and his 
sxit path is hSh2A3A4E2. 

Ths currant asthod for choosing a path is, of courss, not psrfsct. 
Nsascrinr of ttim sffsct of iapsdiasnts by adding othsr pssssngsrs' linss of 
travsl was chossn for its siapLicity and will not, in gsnsral, provids an 
optiaal distribution of pssssngsrs aaong availabls sxits. Ths ration uls 
bshind using ths blocking passsn^ar's tias of aovsasnt as a tias dslay is 
bassd on ^s assui^tion that ths passsngsr undsr considsration will havs to 
wait for ths blocking passsngsr to aovs. This assuaption is aost valid whan 
ths blocking passsngsr is doss to, and slowsr than, the passsngsr under con- 
sideration. If ths blocking passenger is far enough away or fast enough, ths 
other passenger nay never get close enough to hia to havs to wait for hia. On 
the other hand, it can be argued that if a passenger sees another passsngsr in 
his path, he dll view this as presenting a delay. 

More inportantly, at present, the choice of path doss not account for 
any avuxdance of fire or fire-related toxicants. Note that the algorithn does 
not require that be synetric, i.e., that * ^j,i* asans that 

Bovenent towards the fire could be discouraged and aovsasnt a%iay froa ths fire 
encouraged. Future versions of the coaputer progran should include soae con- 
sideration of the tenperature difference between nodes. More analysis is 
required to find the best way to do this, ht present, the aetric is 
tine-based; the concept of distance is considered in terns of the tins the 
passenger believes is required to cover the distance to each unit. While it 
is easy to put a nuaerical factor into the coaputer code to alter d^^j as a 
function of the tenperature in nodes i and j, it is not trivial to detemine 
what realistic values for that nuaerical factor should be. Similarly, a fac- 
tor for confusion as a result of either saoke or panic or both can easily be 
inserted into the computer, but, again, detemxning eui appropriate factor is 
not trivial. 
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Efficient us* of th* axit path algorith* requiras aathoda to dataralna 
how oftan axit oaths should ba updatad. It would ha posslbl* to parfom upda* 
tas after every passanqav aova or other chanqa in th* cabin anvironaant, but 
this would qraatly incraas* aiodal run tia*. At prasant, there are two oacha- 
nisas for drivinq exit path updates. The first is an input paraaetar (sa* 
card type A, Appendix A) which specifies a constant tia* interval between exit 
path updates for all oassengers. Usually, if the interval between updates is 
short, significant alterations of exit path occur for only a saall minority of 
the oassenqers on any given update. This aaans a lot of coaputar tia* is used 
recoaputinq oaths which have not changed. To solve this problea, a second 
method for updating axit paths is provided. This consists of identifying cer- 
tain nodes in the aircraft configuration as decision nodes. When a passenger 
moves into one of these nodus his exit oath is recomputed. Decision nodes are 
places where a passenger has a choice of ways to go, e.g. , row and column 
aisle intersections. (See card type B in Appendix A for inpct descriptions). 
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SECTION 5 
CONCLOSIOH 

The ODRI FTREVAC does not reoresent a finished product. There are two 
bread areas of activity required before the aodel can be depended upon to 
fulfill its stated objective. These areas can be described as oodel valida- 
tion and nodel feature refineaent. The aodel validation work is the nore 
inportant of the two; in fact, it is probable that the pursuit of the aodel 
validation will suggest the direction of aodel refineaents. 

5.1 MODEL VRLIOATICni 

Software testing activities have the dual goals of aodel verification 
and validation. Ne consider verification to be the process ifhereby the ooa- 
puter code is verified to faithfully iapleaent the aatheaatical acdel of the 
siaulation, i.e., where we insure the code is doing what we thought we told it 
to do. vie regard validation as the process whereby we insure that the aodel 
produces an acceptable approxiaation of the real world behaviors it is 
intended to sisiulate. In this fraaeivork, aodel developaent is viewed as a 
building process in which we continually atteapt to iaprove our approxiaation 
of the real world behaviors, expand upon the nuaber and kinds of behaviors 
siaulated, or both. 

The oresent version of FIREVAC has undergone considerable verification 
testing, but very little validation testing. ‘Rie oriaary reason for this is 
ladc of data. Coainsky (see Reference presents a data base resulting froa 
a review of impact survivable oost crash fire accidents. In only a few of 
these is there any data relating to egress rates, and in none ot then is there 
any breakdown of speeds of movement with respects to features of the aircraft 
configuration other than exit chosen. As described in Section 2 , even the 
CAMI tests do not provide adequate data on speed of movenent within the cobin. 
The situation with regard to data on how various toxicants combine to degrade 
passenger movement is even worse. A realistic '^lidation scheme for the 
FIREVAC model must first concentrate on validating the egress simulation in 
the absence of fire, vie must attempt to obtain data for such validation froa 
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manufacturers' certification tests, by analysing videotapes and i4\atever other 
sources of information are available. The ?IREVhC model was purposely 
designed to be heavily dependent «^n input parameters which describe 
passenger movement. With enough data, me should be able to adjust these para- 
meters to obtain a "good” evacuation simulation, where "good* is defined as 
replicating egress rates from emergency evacuation tests. 

The problem of validating tiie poet crash fire scenario is much more 
involved. First, analysis of such accidents does not lend itself to classifi- 
cation. Each accident has so many unique features that a generic class of 
parameter <tescriptions cannot be formulated, i.e., each accident must be 
treated as a special case. Extant descriptions of fire spread, cabin debris, 
passenger conditions, etc. are inadequate. Furthermore, the model's *"iran 
response factor and Fq concentrations are only crude representations of toxi- 
cant effects. Even so, the model can provide a relative measure of hazard for 
different post fire crash scenarios. 

5.2 MODEL FEATORE REFINQfEMTS 

At present, we envision refinements and the inclusion of additional 
capabilities in each of the model's three nodules. Activities under con- 
sideration include: 

CEH 

1 ) Smoke could be included as a function of both time and cabin 
position. 

2) Fire scenario input needs refining. This could include the ana- 
lysis of data from the HASA ilouston fire tests (see Reference 5) 
to refine the approximation of cabin atmosphere data for all 
cabin nodes from the test data, and exploring the possibility of 
using more sophisticated fire ondels to produce input, we should 
note that the techniques of sophisticated fire models (e.g. , 

Metre Dame's DNDSAFE) which are PDE (partial differential 
equation) solvers, require amounts of computer time and space 
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which preclude trying to incorporate those techniques in our 
Bodel. He see the deyelopeent of sets of "representative 
situations* frosi test data and/or PDE aodels as our best alter- 
native. 

3) Inclusion of crash-related cabin debris could be included. 


HFM 

1) The calculation outlined in Section 3.1 could be augaented 
considering each paissenger's body mass and respiration rate. 
This would result in: 


FoCt) 




, T_ 

0 I- Q 


3.61 


R C. 

. -I 1 
^ d. 


dt 


where F^(t) is the fractional incapacitation dose at time t, (-) 


is the airisient gas toqterature (C) 


3.61 


H 

R 

V 

d. 


is an eapirical constant (Crane) (c*'"' -sec) 

is the aabient concentration of the ith toxic gas (ppa) 

is the passenger response factor (-) 

is the body mass in ga 

is the respiration rate in al/sec 

is the incapacitation dose in PPM al/ga. 


2) The F^ calculation could be replaced with the concept of a short 
tern lethal Unit, passenger incapacitation would be assuned 
instantaneous upon absorption of a specified lethal dose of any 
toxicant. 


3) Replace an additive with the aaximua fractional incapacitation 
dose of the indi/idual toxicants as absorbed at time t. 


PEN 

1 ) we need to obtain and evaluate data from emergency evacuation 

tests and use the results to better define passenger speed para- 
meters . 
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2) >ie have to inprove the dioica of exit path process . nils 
includes alterinq the distance function to reflect the presence 
of firs and seoke as %iell as incorporating soee sort of confusion 
factor due to panic. Also required is refineaent of the exit 
path update criteria; %dien should a passenger dianqe his aand 
about choice of exit? 

3) Vie need to consider the inposition of delays caused by panic or 
confusion. 

4) The possibility of a <dvange in contender priority logic (as 
discussed in Section 4.1) naybe desirable. 

Model inorovenents such as those listed above should be given priority 
as a function of their oromise for support towards the model objective. The 
model objectiw itself should be refined to detemine how the nodel is to be 
used, and vdiat the specific purposes of exercising the nodel are. 
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The FIREVAC nod«l requires several different types of input data and 
produces a variety of outputs. In the interests of modular structure the coai<» 
puter program has been designed to read from, and write to, different logical 
units. Each logical unit is associated with a VAX 11/780 file, and each VAX' 
file can be considered to contain a distinct type of input or output data. 
Figure 7 shows a s^UBple VAX command file, a list of the steps required to run 
FIREVAC on a VAX. The procedure illustrated in Figure 7 assumes that am 
executable version of the FORTRAN code is in memory (in this case chat version 
is called FIRESND) . It further amsumes chat the required input files are 
available. The input files aure assigned to the FORTRAN input device numbers 
11 >16 (FOR011 - FOR016). The specific variable formats for each file are 
included below. Briefly, the input files are: 

1) A control file, assigned to FOROll. 

The control file contains model control parauneters such as model 
cycle time updates. See Table l below. 

2) An aircraft configuration file, assigned to FOR012. 

This file contains the physical description of the aiircraft with 
definitions of each node (e.g., seat, aisle, etc.). See Table 2 
below. 

3) A passenger position file, assigned to FOR013. 

Ihis file defines Che number of passengers and their initial loca- 
tions. See Table 3 below. 

4) An aircraft atmosphere file, assigned to FOR014. 

This file contains the cabin atmosphere profile. This file is 
required even if the simulation is to be run without an atmosphere. 
The model requires the number of toxicants and toxicant sampling 
points - which would be set to 0 in the no atmosphere case. See 
Table 4 below. 
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S! Input Fil«s: 

•ASSIGN DATA: CTRL180A. OAT FOROll 
•ASSIGN DATA: ACISOCMCL. OAT F0S012 
•ASSIGN DATA: PAS8180. OAT F0S013 
•ASSIGN DATA: H14A. OAT F0R014 
•ASSIGN GNOVe: FHYSINDAT. OAT F(3R019 
•ASSIGN DATA: GOEV. DAT F0R016 
•! Output Fil««: 

•ASSIGN FATHI80. OUT F0R021 
•ASSIGN SNAPiaO. OUT F0R022 
•ASSIGN ATNQSISO. OUT F0R023 
•ASSIGN PUOTiaO. OUT F0R024 
•ASSIGN SUNRYISO. OUT F0R0S9 
•SUN FIRE9NO 


! 0C9 — SERIES 180 cat*. 
! FIRE9ND* atamtphtrt. 


Figure 7. VAX Command File to Run the FIREVAC Simulation. 
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TABLE 1 

CONTROL FILE 





FOROin •• 

Paaaangtr •frcaa alMulatlan 
atark klaa. Uaa4 ka tnlklallia 
TOFi-H : Paaaanfar tiaia oF laak 
•avaaank. 

Blaulaklan and ti»a. Haiiaiua 
allowabla clock kl*a. 

Clack u^Faka Incraaank. Uawally 
oak ko 10 (.01 aac. a) aa kfclo la 
an orbar of aa|nlku4a aaallar 
kkan klia Faakaak ^aaaanfor aeva- 
•ank klaaa. 

Tlaa Inkarval bakMoan oukpwk 
anapabaka. 

Tlaa Inkarval bakMoan aknaa^bara 
w^baka calla. Tba alrcrafk 
aknaoFbara la aaawnab canakank 
bakataan upbakaa. Bhawlb ba aak 
ka a valwa graakar kban kba 
aagackab alaulaklon karalnaklon 
IF no Flra-ralakab aknoagbarlc 
oFFacfca ara balng alaulakab. 

Tlaa Inkarval bakwaan gaaaaniar 

canblklan (HFHl wgbakaa. Bhawlb alaa 
ba aak ka a valwa graakar kban 
kba aigackab alnwlaklon 
karalnaklon klaa Far kba no- 
Flra (no abnotFboral acanarla. 




TABLE 1 
CONTINUED 


A-7 

JPSUBO 

lnt>q«r 

1-7 

16 


Dabtifflng g«ra«atar. UaaB to 

grlnt akaioaRharlc affacta an tha 
JPBDBOtli gaaaangar. .If no awch 
Butgut la Jaalratfi alioulB ba aat 
to 0. 

A-B 

I X 1 TUP 

Intagar 

1-6 

16 

•till aac 

Tina Intarval batiaaan goaaanior 
ailt gatb ugBatoa Marnlnf — 

tbla garanatar baa a larga 

affact on grograai run tlaio. 

A-9 

REACT 

Raal 

1-6 

E6 O 

aaca. 

Initial oaaaangar raoctlon tina 
Tbla volua la 4»§rm4m4 for aacb 
gaaaantor aa cabin conbltlona 
dlctata. 

A- 10 

lOaCRB 

Char 

1-B3 

73A1 

- 

Run Baacrlftlon of currant 
alnulatlon. 
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TABLE 2 

AIRCRAFT CONFIGURATION FILE 


CJ 


Vmr labia 


lacord Naaa 

Tgpa 

Col a. 

Foraak 

Uni to 

Daacr Ipklon 


aa 

READ IN 

SUBROUTINE 

AC IN Ifroa 

FOROI2) oa 

B-1 lACrVP 

Char 

1 -to 

lOAl 

- 

AIrcrafk tgpe. 

B-2 MUIDE 

Intaper 

1-3 

13 

- 

Length of aircraft 1 In 
nuahar of no4aa> 

NLUNO 

Intagar 

4-6 

13 


Nldkh of aircraft lln 
nuahar of nodaa). 

B-3 JCNOX 

Integer 

1-3 

13 

- 

Coluan nuahara in which 


•l«l«b ar* locatarf 

Laafc antry “ O. 


IRNOX Intager 1-3 13 Row nwabara In which 

alblta art locakai. 
Laak anirg “ O. 

B-9 <Una kgpa B-9 card ahoul4 ha raad In tor aach ot kha MM IDE ■ MLOMO 

alrcrafk na<laa l a.. NR ahould varg froM 1 ko NLONO aa NC varlat 
froa I ko MMIOE ) 



Row nuahar of currank no4a. 

Coluan nuahar of currank no4a. 
Currank noda'a agaclftc notfa kgpa: 

1 “ row alala 

2 “ col alala 

3 - aaak 

4 “ ailk. kgpa A 

9 ailk. kgga A OM 

6 “ aliga. kgpa A 

7 - alMa. kgpa A OW 
B ak In/ai kar lor 

NTVPE goaa inko NDETVPINR. »4C I 
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TABLE 2 
CC»iTINUED 


LJ 

OV 



Nce 

Integer 

21-23 

12 


B~6 

NUrINUE 

Integer 

6-a 

13 

- 

B-7 

XBARi IHi 

Real 

1-7 

F7. 4 

- 


a< IN> 

Reel 

9-19 

F7 4 



TrtAK4 INI 

Reel 

17-23 

F7 4 



iniNt INI 

Reel 

29-31 

F7. 4 

- 

B'8 

NDECPT 

Integer 

1-3 

13 

- 

B-9 

<A B 
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le read 

for each 

doc telon goint 


lOECPI 
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1-3 

13 

- 


<MPECPI 1 






juetpi 

Integer 

4 6 

13 
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NC8 goa* into NCABEiIM.NC) 
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UN«1>NUHN0E> 
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travel Ing tbrougb noO. 
<iH«i.HuriNDe> 

Saagla »tan4ar4 bavtatlan 
for no4a. 

(IN-t.NUHNOrt 

HadaiMO •goail througb na4a 
4 IN-I.NUMNDEI 

ninu*u« sgaa4 ttirougb no4a 
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INOECPTI) 

Row location of 4ac iolon 
goint 

Cotuan location of boclolon 
foint. 
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CONTINUED 


Vai-lakI* 
Rvcord Nsai* 


Ty^> Cola. Foraiat Units 


Ooscrlftlon 


C-l 
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• • 

Intagar 

RfcAO IN «• 
1 -a 

*BROUTINE 

la 

EXITIN 

(fi-oai POROiai as 

Nuaher of salts In the 

c-a 

(Una 

typa C-a 

card should ha read 

in tor 

aircraft. 

each exit (1-t.NEXITSM 


LOCRXI 1 > 

Intagar 

1-3 

13 

_ 

Row location of aslt. 


LOCCXdl 

Intagar 


13 

- 

Coluain location of oatt. 


IXIIPPI 1 » 

Intagar 

7-9 

13 


Note that this value should 
either HROTXT or LEFTXT. 
Passonier forceftlon of 


TOPEN( 1 ) 

Real 

la-ta 

F7. 3 


salt status: 0 “ ogan> 

1 closed. 

Tlaia It tatas for salt to 

ha ofenad. 
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TABIJi 3 

PASSBNGER POSITION FILE 


UJ 

00 


N»c or 


D-1 


0-2 


0-3 


Vorlakl* 


Naae 

Igpa Col ■. 

Foraut 

Unit! 


•• REAO IN 

BURROUTINE 

PASS IN 

1 AC 1 VP 

Char 1-10 

lOAl 

- 


NUtIPAB 

Intvtar 

1-3 

13 


IPASa 

Intagar 

1-3 

13 


NR 

Intrgar 

9-7 

13 


NC 

Intrgar 

9-11 

13 

- 


Dotcrl^tion 


Ifroa F(mOi3l •• 


Aircraft tyf* 

Not* that hhli tliawltf aatch 
racar4 l-l af tha aircraft 
Infaraatlan flla 

Nwabav- af aaatf’f*'** •" 

tha aircraft far thla 
•iawlatlan 

Paatangar nuahar 

Initial ra« lacatlan af 

gatoangar 

Initial caluan lacatlon of 

aaaaangar. 

tlata that thaaa valwat ahouia 
arahahig inficata a aaat lacatlan 
at Patcriht# on cart tgga •■9 

1 a. NED1VPlNR.NO - 3 


9S 
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TABLE 4 

AIRCKAET ATMOSPHERE FILE 


vC 


V#i»r* le 
flacortf Na«* 


E-l 


6-2 


6-3 


6-4 


6-6 


T»6* 


Cu 1 • 


Fomat 


Units 


OascrlFkton 


NO ABES 


aa read in 8UIR0UIINE ArHSIN Ifroa FfMOMI 
Intagar 1-4 14 


<6 2 aii4 E-3 
NANOAaC I. 10) Cliar 


OABIOCI ie> 


NTEHP1 


Raa I 
Intagar 


<Aii 6-3 card li 
6-3 NOASFtI) Intagar 


l-IO lOAi 


1-4 


FI3 6 
14 


3-64 


14 


INANE Cliar r-20 20AI 

atffarant toilcanta In tha 


PPN»aac 


Nu«kir at gaitt to 

Og thta alMulation 


Naaia of aacO gaa. 

Not# that orOar oF gaaaa 
•gaclFla4 aiuat raaaln canalatant 
with tha orFar uaaO ta Ingut gaa 
cancantrablan valuaa halo«. 

I l-l. io> 

Incagacltatlon Oaaa conatanta 
IPPnaKC) For aach gaa 

NuaAor eF taaigaratwra aaaigling 
galnta. 


Nuaihar oF aaaigltng gotnta 
aaaaclata4 with aach 
iiiOiviaual gaa. 

Daacrigtion aF Fata coaMoant 
UaaF to Faliailt valwaa Far 

ataaoghara graFlla Input Flla 

Tina caorFInata Far all alF 
taaiparatura anO gaa eoncfn- 
tratlan valuaa 
Bhauli ha <- BIHSI 


«S 


E-7 


1 1NOI II 


Rtal 


1-iO FIDO 
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TABLE 4 
CONTINUED 


E-a 

E-9 


E-IO 

E-ll 


4E-8 thru E-IO cards ara ra^ulrad only id NTEHPf>0> 


INAME 

CHAR 

1-30 

80AI 


Ooscriptlon of data coMMont 
usad to dlstlnpuloh tonporaturo 
valuot In tho data fllo 

(Card 

typs E-9 

Is ropoatod until 

all toMporatura i 

lanpllni points 

mrm 

Input. (1 

-l.NTEMPTl 

Data 

for 4 data points 


oach 

card — 1 . 0 

1 fornat 

Is 4(14 

. I4.F10 01) 


iTtrTdi 

Intofor 

1-d 

14 

- 

Row coordlnato of oach 



19-33 

37-dO 

98-90 




jimiii 

Intofor 

8-a 

14 

- 

Colunn coordinata of aach 



33-86 

41-44 

99-43 




THPOLDI 1 > 
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9-ia 

no. 0 

dot s C 
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49-94 
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IFTth sanplint point. 
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1-10 

FIO. 0 
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taadaratura valuat iunit lanftb 
« nada ltn|lh ar widthl 
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•asas in tha data Mis. This 
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of tha |as. Mata that orasr In 
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CONTINUED 


E -la 


(Card tiira E 18 It rt^takad until all taapllnf for IQtt» t** 

doacrikad Data for 4 data folntt It attwaad on oack card--- 

1 a . fornut Lt 114. 14. f 10 0> > 


lOABPI 
( IIM. 

lU) 

IntoMOf 

1 -4 

14 


JOABP 1 
( Iff. 

101 

Intofor 

S B 

14 

- 

0A8UI 0 
< IPI. 

IU> 

Naal 

9-lB 

FlO 0 
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Rom coordlnata of IRtk canconkratlan 
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taotfllnf folnt for lOtk fat. 
<IPT-i,NOABIO(iei* lO-I.NOABBB) 
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Intorf olaklanl canconkratlan 
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(NUIE. Old taaioaraturo and fat canconkratlan vaiuot aro road In taltk 
location of taofllrif foinkt Nom valuoo aro ooauaad to fco at taaa 
taapllna rnlnl* In taaia ardor. ) 
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1-10 flO O 


nodoo 
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TABLE 4 
CONTINUED 


V«rlakl» 
Record NaMo 


Typo Cult For«at Unit* 


Datcr lotion 


r> 


F-l 

F a 


F-3 

F-4 


READ IN SUBROUTINE ATM8R0 tfroM FDROI4I #• 


IF-'l thru F-6 type card* are repaatad For aacte tlaa I ta«para%urat 
pat concantrat ionl lahara kha latk tliao thould ba iraakar than or 
at|i>al to anticipatad run taraiinakion tti<..l 


INANE Char lao 30AI 


Oatcrlpklon oF data cedNoant. 


riNNEN Real 


l-IO FIO O aaca. 


Tl«a caardlnata Far all 'naw' 
ta*paratura and t** roncantratlon 
valuat. 


(F-3 and F~4 t«tpa rardt ara rapulrad ontp IF NTEMFT > 0. I 


INANE 


Char 


1-ao 


SOAl 


Oatcrlpklon oF data catMaonk. 


TNFNEW(I) Real 


l-IO BFIO 0 da|. C 

1 1 - ao 

UI-3U 
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41 30 
31 -AO 
61-/0 
71 no 


ll-l.NTEHPTI 
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aatipllnf point. Ordar oP valwat 
abould ba kba tana at In TTrtN.D 
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For each pat For Mhich N0A8PTU0) > 0. Tha F-B card abowld bo 
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Fur all tawpllnp polntn I 


•n 5 

S z 

€2 
> o 
^ m 

3a 


F - 3 


INANE 


Char 


l-:*0 OOAI 


Oatcrlpklon oF data coaMuant. 
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CONTINUED 


F-6 OAHNtW 

( ir r. iti> 


1-10 8FI0 O PP«« 

1 1 ao 
at- :io 
31 40 
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N«M vAlu* for lOtfc ••• coneon- 
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5) A passenger descripcion file assigned to FOROIS* 

AC present this file is noc used, it is intended to provide 
detailed passenger descripcj.ons (e.g., sex, age, body mass, respira- 
tion rate) as model features requiring such data are inplenented. 

6) A file of Gaussian deviates assigned to POR016. 

This file provides a normal distribution of passenger charac- 
teristics (speed) for the present model implementation. See Table 5 
below. 

itodel output files are: 

1) Passenger movement details, assigned to F0R021 . 

This file is intended primarily for debugging purposes. It lists 
every passenger move on a node-to-node basis. Debugging switches in 
the model can be set to greatly expand this file, giving exit path 
information or examining the results of individual subroutine 
calls . 

2) Snapshot output; assigned to FCR022. 

This IS the graphical representation of the aircraft interior and 
passenger positions (see Figxires 3 to 1 6 of Appendix B) 

3) Atmosphere output, assigned to FOR023. 

This IS a record of the cabin's interior atmosphere. Here again, the 
volume (and amount of detail) can be controlled by debugging 
swi tches . 

4) time plot data, assigned to FOR024. 

Provides a sec of data points with time as the abscissa, and the 
fractional incapacitation dose of a selected passenger as ordinate. 
This output can be used to generate a plot of vs time. 

5) Simulation summary output, assigned to FOR015. 

This gives the time of evacuation for each passenger, the number of 
passengers evacuated through each exit and last passenger's time out 
for each exit. See Figure ’ ^ in Appendix 3. 
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TABLE 5 

GAUSSIAN DEVIATES FILES 


Var 1 ab I e 

Record Nam» T^pu Cole Foraab Unit* UascrtFtlon 

a> NFAU IN 8UBROUIINE INUSFC Ifroa FCM016I •• 

<0 I Cards are rapaatad until NOEV Oausslan davlalaa liava been 
read In At praaank NOEV " 350 I 

O-l 0 REAL 1-a 9(FI3. 01 Oausslan davlakas 

13-34 
39 36 
37 40 
4V-60 


16. 

01 
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APPENDIX B 
SAMPLE OUTPUT 


ORIGINAL PAGE IS 
OF POOR 


This appendix contains a comparison of a no fire scenario with a fire 
scenario. The airplane consider^ was a OC-9 5180. For the fire scenario the 
atmosphere data used was taken from test case 24 of reference 5. This output 
is not intended as a prediction of actual evacuation, but is provided solely 
to illustrate typical model output. Figures 8 to 16 show passenger placement 
for both scenarios at 20 second intervals. Figure 17 shows a comparison the 
summary outputs for both cases. The exit numbering scheme used is shown in 
Figure 8. 
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Figure 8. 

Initial passenger placement 
in both Fire and No Fire 
Scenarios. 
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Exit #4 
Exit #6 


Exit #10 
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Figure 9. 

Passenger placement at 
t = 20 seconds 
No Fire Scenario. 
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Figure 10. 

Passenger placement at 
t = 20 seconds 
Fire Scenario. 
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Figure 11. 

Passenger placement at 
t » 40 seconds 
No Fire Scenario. 
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AT rtnt« 40. 00 
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Figure 13. 

Passenger plact-ent at 
t = seconds 
No Fire Scenario. 
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Figure 14. 

Passenger placement at 
t =60 seconds 
^ire Scenario. 
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Figurt 15. 

Passenger placement at 
t = 80 seconds 
No Fire Scenario. 
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aciT 

Na OUT 

TOTAL TIME 

1 

48 

63.603 SEC. S 

2 

32 

63. 462 SEC. S 

7 

83 

99.869 SEC. S 

10 

IS 

30. 124 SEC. S 


SIMULATION START* 0 MILLISECONDS DELTA T* SO MILLISECONDS 
ACCELERATION FACTOR* 0. 290 
OC-F 3180 

OCY S-180 FIRE9N0 COMPLETE ATMOS (IATMDT*I0TSP0*9 SEC. S) 

(a) Fire Scenario (Atmosphere and Human Factor Updates Every 5 seconds). 


EXIT NO. OUT TOTAL TIME 


1 

46 

63. 823 

SEC. S 

2 

39 

64. 023 

SEC. S 

7 

76 

83. 493 

SEC. S 

10 

17 

64. 140 

SEC. S 


SIMULATION START* 0 MILLISECONDS DELTA T* SO MILLISECONDS 
ACCELERATION FACTOR* 0. 290 
OC-9 S180 

DC9 S-180 FIRES Mo Ataos. (IATHDT-9'*/IDTSPD-9's/lXITUP*l00000) 

(b) No Fire Scenario. 


Figure 17 


Comparison of Exit Times Between a Fire amd No- fire 
Scenario. 
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model SUBROOTINES 
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This appendix provides the details of the eodel's subroutine structure. 
The model is first divided into its specific procesr:es and then each process 
is in turn broken down into its components until the model is described in 
terms of its individual FORTRAN subroutines. 
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FIREVAC VERSION 5 PROCESSES 


Atmos phera Update 

Exit Path Update 

Clock and Configuration Update 

Human Reaction Factor Update 

Input/ Ini tialization 

Passenger Movement 
Snapshot 


Calculates cabin atmosphere from input 
atmosphere profile. 

Determines "optimal* exit path for a given 
passenger. 

Keeps tunning sumulation clock controls. 
CQ( and HFM updates. 

Calculates human reaction factor and 
adjusts passenger speeds. 

Reeds in required aircraft and passenger 
data. Sets initial values for the model's 
dynamic variables. 

Simulates passenger movement from seated 
position to exit. 

Produces graphical output showing cabin 
interior and passenger position. 
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FIREVAC VERSION 5 SUBROUTINES 


ACIN 

ADJCNT 

ATOSIN 

ATMSUP 

ATMSRD 

CLKUP 

CNFGUP 

CNTRIN 

CONTND 

EXITIN 

EXTSIM 

INDSPC 

INTRAFA 

ISTNBL 


Aircraft input routine reads aircraft configuration data, type, 
width, length aisle locations, etc. 

Finds all nodes vdiich are adjacent to a given node and which can 
be occupied by a passenger. 

Atmosphere input/initialization routine. Reads in profile of 
atmospheric conditions and initializes other atmosphere related 
variables. . 

Updates the atmospheric conditions inside the aircraft for a given 
time. Used first to set initial cabin atmosphere parameters then 
called during the sumulation to update them. Calculates toxicant 
concentrations to reach node in the aircraft cabin. . 

Reads one "new" set of toxicant concentrations from atmosphere 
input profile. 

Clock update routine. Updates simulation clock. 

Configuration update routine. Updates the condition of the 
aircraft. 

Simulation control data input routine: start time, update times, 

etc. 

Finds all possible contenders (passengers wanting to move into) a 
given node. 

Exit information input routines, reads number of exits, locations 
whether exit is open or closed, and time to open exit, etc. 

Controls cycle between PEM, CEM, auid HFM. 

Individual Specification. Reads a set of Gaussian deviates and 
uses them to calculate each passenger's speed through each node 
type. 

Linearly interpolates between "old and "new" toxicant values from 
the atmosphere input profile. The interpolation is on time depen- 
dent sampling point temperatures and gas concentrations. Spatial 
approximations are handled by subroutine WGTAVE. 

Is target node blocked? Routine decides whether a passenger's 
target node (next in exit path) is blocked. 


61 



METRIC 

MOVE 

NUMXIT 

PASSIM 

PHYSIN 

PPOSUP 

PRTYMV 

SNAP 

SPDCHG 

SPEED 

SUMRY 

TRCRTE 

TTOMOV 

WGTAVE 

XITPATH 
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FIREVAC VERSION 5 SUBROUTINES (Continued) 

Measures "distance" between tvio nodes, for use in exit path 
calculaticps. 

Move routine moves the passenger into target - ode. 

Service function which returns a value of zero if a given node is 
not an exit, and otherwise returns the exit number. 

Passenger input routine reads initial (seat) passengers locations. 

Initialize passenger physical characteristics (at present only 
sets » 0 for all passengers). 

Passenger position update passengers who are able to move at pre- 
sent clock time will be advanced along their exit paths. 

Priority Move. Determines which of the contenders for a given 
node has priority. Present criterion is longest waiting time. 

Snapshot routine generates a rough presentation of the columns and 
rows to give a snapshot of where each passenger is located in the 
aircraft at a given time. 

Uses human reaction factor to adjust passenger speed through node 
types for a given passenger. 

Updates the fractional effective dose, and then calculates a speed 
factor for that passenger. 

Summary routine prints summary of exit output data. 

Traces the path to the nearest exit for a given passenger. 

Time to move? Routine decides whether enough time has elapsed 
since last passenger move, to allow passenger to move again. 

Uses sample point data of the input atmosphere profile to calcu- 
late a weighted average value for a given toxicant and a given 
node. 

Finds optimal paths to open exits for a given passenger. 
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